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Abstract 

Different momentum space Faddeev-like equations and their solutions for 
the radiative pd-capture and the three-nucleon photodisintegration of '^He 
are presented. Apphcations are based on the AV18 nucleon-nucleon and the 
Urbana IX three nucleon forces. Meson exchange currents are included using 
the Siegert theorem. A very good agreement has been found in all cases 
indicating the reliability of the used numerical methods. Predictions for cross 
sections and polarization observables in the pd-capture and the complete three 
nucleon breakup of ^He at different incoming deuteron/photon energies are 
presented. 
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I. INTRODUCTION 



In the case of few nucleoli systems it is nowadays possible to compare precise experimental 
data and theoretically well controlled predictions. For low energy processes with three 
nucleons (3N), theoretical results can be obtained for any given realistic nuclear interaction. 
This makes such systems an important tool for the investigation of the nuclear Hamiltonian. 
In momentum space the formalism of Faddeev equations has been used to obtain bound and 
scattering nuclear states [1]. A very good agreement between theoretical predictions and 
experimental data was obtained (e.g. [2]). An equivalent description was also obtained using 
configuration space variational methods [3] . After such encouraging results in pure nucleonic 
systems had become available also weak and electromagnetic processes with three nucleons 
were studied in the same scheme. In a series of papers we showed the results for electron 
induced processes [4] , proton-deuteron radiative capture [5] , muon capture [6] and 3N bound 
states photodisintegration [7,8]. It was found that all dynamical ingredients are important: 
final state interactions play a significant role and the clear effect of 3N forces is noticed. On 
top of that also the addition of meson exchange currents changes predictions in a significant 
way, making the analysis based only on the single nucleon current mostly meaningless. 
Inclusion of all those components allows for precise predictions on different observables, like 
cross sections or asymmetries and consequently the door is open for investigations of other 
physical issues e.g. neutron electromagnetic formf actors [9]. 

For very low energies hyperspherical harmonic expansion methods were used by the Pisa 
group [10] and a nice agreement with our results was observed. The total photodisintegra- 
tion cross sections were also calculated by the Trento group [11] using the Lorentz integral 
transform method. We compared our results in a common paper [12]. Also the Hanover 
group presented results on photo- [13] and electro- [14] disintegrations. In their approach the 
A degree of freedom and corresponding nuclear currents are taken explicitly into account. 
These predictions are also in qualitative agreement with ours. 

In the Faddeev scheme different formulations of the three-body problem are possible. 
In this work we would like to compare different ways of obtaining transition amplitudes 
for two- and three- body photodisintegration of 3N bound states. The comparison between 
predictions based on the different formulations is a strong test of the used numerical methods. 
The numerical calculation of the three body continuum is a non-trivial problem and the 
possibility of testing different formulations deserves thorough investigations. Up to now even 
without 3NF's, no direct systematic comparison of exclusive three-body photodisintegration 
cross sections between different approaches was performed. Including 3NF's the situation is 
even worse. To the best of our knowledge, no other collaboration has done calculations with 
explicit 3NF's. Therefore, an internal comparison is of utmost importance. It also provides 
useful information on the efficiency of the different formulations for practical calculations. 

The calculations presented in this paper are based on the AV18 NN potential [15] alone, 
and combined with the Urbana IX 3N force [16]. Since the scope of this paper is not 
an investigation of details of the electromagnetic current operator, the same model of the 
current is used in all investigated formulations. The single nucleon current is supplemented 
by some exchange currents included by the Siegert theorem. This approach is described 
in [5] in more detail. 

In Section II we describe three ways of obtaining the transition amplitude for the radiative 
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Nd-capturc (and equivalently for the two-body photodisintegration of the 3N bound state) 
and two methods for the three-body photodisintegration of the 3N bound state. In Section 
III we compare predictions based on those methods. We summarize in Section IV. 

II. THEORETICAL FRAMEWORK 

In this section we would like to describe three different methods used to generate transi- 
tion amplitudes for the Nd-capture and the two-body photodisintegration of the 3N bound 
state. We also show how to build the transition amplitudes for the three-body photodisin- 
tegration. 

The nuclear matrix element for the two-body photodisintegration of the 3N bound state 
I *6) is 

^ (*Sv"j I J, I *^>) , (1) 

where (^'^j | is the final scattering state. In the Faddeev scheme it can be presented in the 
form [7] 

TV;-^ = (V^i I (1 + P)j, I , (2) 

where {ipi \ is a Faddeev component of | ^'^j) and is the electromagnetic current operator. 
P is a permutation operator defined as a sum of cyclical and anti-cyclical permutations of 
three particles 

P = -P12-P23 + -P13-P32, (3) 

where Pij interchanges the i-th and j-th nucleons. 

The Faddeev amplitude (■^i | obeys the Faddeev-like equation [17] 

(^1 I = I +(^^1 I [PtiGo + (1 + P)vl'^GoihGo + 1)] 

= ((/.I I +(^^1 \K, (4) 

where | (pi) is a product of the deuteron state and a momentum eigenstate of the spectator 
nucleon. Further, v}^^ , Gq and ti are a part of the 3NF symmetrical under exchanges of 
nucleons 2 and 3, the free three- nucleon propagator and the two-body t-operator acting in 
the 2-3 subspace, respectively. Thus 

(V^i 1= (01 I (1 - K)-' (5) 

and 

N;^'^{<I>i\{1-K)-\1 + P)j,\^,). (6) 
Defining the auxihary state | U) 

|^)^(l-K)-i(l + P)j^|*,) (7) 
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one gets 



= ('/'I I U) . (8) 



According to the definition (7) the state | U) fulfills 

I [/) = (1 + P)j^ \^,) + K\U) . (9) 

Inserting K this reads 

I [/) = (! + P)j^ I 

+ [PhGo + (1 + P)vPGo{hGo + 1)] I U) . (10) 

This form of the kernel with P standing to the left causes unnecessary complications since 
the deuteron pole in ti appears as smeared-out into a logarithmic singularity [1]. We avoid 
that by reformulation of Eq.( 10), as is shown below. 



A. Methods 1nn+3NF and I^vat 
Denoting | x) = (1 + -P)Jm I ^^'^ introducing the auxiliary states | U') and | U"): 

I U') = hGo I U) (11) 
I U") = V^GoihGo + l)\U) (12) 

one gets 

+ (01 I (1 + P) I U") . (13) 
The states | U') and | U") fulfill the set of coupled equations: 

I U') = ^1^0(1 + P)j^ I ^b) + hGoP I U') 

+ hGo{l + P) I [/") (14) 
I U") = vl'^Goil + hGo) I (1 + P)i^ I 

+ vl^^Go{l + tiGo)P\U') 

+ v}'^Goil + hGo)il + P) I U") . (15) 

Solving numerically the set of Eqs (14)-(15) and using Eq. (13) one gets the transition 
amplitude N^'^. In the following, the results based on the Eqs (13)- (15) will be denoted as 
"method 1nn+3Nf" 

In the case when only the NN interaction is used {v}^^ = =^| U") = 0) Eq. (14) 
simplifies to 

I U') = hGo{l + P)j^ I ^b) + tiGoP I U') (16) 

and 
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= (-^1 I X) + (01 I P I U') (17) 
This will be called "method Inn" 

In all our methods the inhomogeneous integral equations will always be solved by iter- 
ation and consecutive Pade summation. In the numerical implementation it is important 
that in both methods, during the iterations of the set of Eqs. (14)- (15) or Eq. (16) the 
permutation operators from the integral kernels act only onto the t— operator or the 3NF's 
matrix elements, which stand at the very left in the driving terms of Eqs.( 14)-(16). We 
work in a partial wave decomposition and the presence of the nuclear interactions (in ti and 
V4) enforces that only channels with relatively small partial waves are important. Thus the 
P operator which acts upon ti or V4 can also be taken using a relatively small number of 
partial waves. 

Having solved the set of Eqs. (14)-(15) one can also obtain the amplitude for three-body 
photodisintegration [7] 

= (03iv I (1 + P)j^ I + {hN I (1 + P){\ U')+ I U")} . (18) 

Here (03iv | is a product of momentum eigenstate describing three free nucleons and anti- 
symmetrized in the 23 subsystem. 

B. Methods 2jvAr+3ArF and 2nn 

The second method, which we will denote as "2]sfN+3NF" has been presented in detail 
in [7], where also some predictions for the pd-capture and the two body photodisintegration 
were shown. However, for the purpose of completeness we briefly describe also this method. 
Using the identity 

1 + P = ^P(1 + P) (19) 
and introducing the auxiliary state | U): 

I U) = tGo I U) + ^(1 + P)vj'''Go{tGo + l)\U) (20) 
one gets from Eq. (10) 

\U) = il + P)j^\^^) + P\U) . (21) 
Then the Faddeev-like equation for the state | U) is 

I U) = (tGo + ^(1 + P)vi'^Go{tGo + 1)) 
* {1 + P)j^ I *b) 

+ (tGoP + ^(1 + P)v}'^Go{tGo + 1)P I U), (22) 
and the transition amplitudes are 
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AT;^^ = (01 I (1 + P) I I + (<^i I P I C7) (23) 

and 

^r=(</>3ivi(i+p)ij^b) 

+ I tGoil + P)j^ I ^b) 

+ ihN \P\U) + {<P3N I tGoP 

In the case when only the NN interaction is used {V}^^ 

\U)^\x)+PtiGo\U) . 

Using two consequences of the identity (19): 
\ x) — ^P \ x) 2 = P{P — 1) Eq. (25) can be rewritten as 

\U) = ^P[\x)+hGoP{P-l)\U)]. (26) 

This is equivalent to 

I U) = \P I U'") , (27) 

where | U'") fulfills 

I U'") =1 x) + hGoP I [/'") . (28) 

This equivalence can be easily seen by iterating Eq. (26) and Eq. (28). Therefore one obtains 
from Eq. (8) and Eq. (27) 

= \ P \ U'") ■ (29) 
The transition amplitude for the three-body photodisintegration is given by [7] 

= ^(<^3^ I P I U'") + lihN I tGoP I U'") . (30) 

The results based on Eqs. (29) and (30) will be denoted as "method 2jvjv." In both meth- 
ods 2is[N+3NF and 2jviv one meets the action of the permutation operator P onto matrix 
elements of the P-operator from the previous iteration. In the actual numerical implemen- 
tation the identity (19) is fulfilled only approximately. Since we work in a partial wave 
decomposition, wc can take into account only a finite number of partial waves. In order to 
fulfill the identity (19) a high number of partial waves has to be used. However, also the 
number of necessary two-body channels increases with the value of total momentum of the 
two-body subsystem. In consequence, this requires a large amount of memory and comput- 
ing time. We have found that the worst convergence occurs in the case of method 2jvAr and 
one needs to use two-body channels with j„tax = 5 during the iteration of Eq.( 28). It is 
due to the lack of the t-matrix in the driving term of this equation. In all other methods 
the t-matrix, which is most active in the lower channels, reduces the infiuence of higher 
two-body channels. Therefore even using the identity (19) one can restrict the number of 
partial waves. 

For all the methods described above, one obtains the transition amplitude A*"™*^ for the 
radiative capture using time reversal in the two-body photodisintegration amplitude Nj^'^. 



U). (24) 
= 0) one gets from Eq. (10) 

(25) 



6 



C. Methods Snn+3NF and Sjvat 



Another possibility to obtain the transition amphtude NJ^"''^ is to calculate it directly for 
the radiative capture process. Then methods used for the elastic Nd scattering can be used 
where the initial Nd channel state occurs in the driving term of the corresponding equation 
[5]. This makes a clear difference to the previous methods. Then the matrix elements of the 
nuclear current comes in the last stage of the calculations after the solution of the Faddeev 
equation is obtained. 

In this case we calculate directly the transition amplitude for the radiative Nd-capture 

n;^ = 1 j, I *Sv1) • (31) 

The Faddeev component | ipi) forming | '^^^}) as 

I *Sv1) = (1 + I V^i) (32) 

is given via 

I V^i) =1 0i) + Gof I 0i) , (33) 

where T \ (pi) obeys 

T I 0i) = I 01) + (1 + tGo)v}'\l + P) I <^i) 
+ tPG^f I (Pi) 

+ {l + tG^)vi'\l + P)Gof I 0i) (34) 

Solving Eq. (34) one gets the amplitude T \ (pi). The next step in the numerical imple- 
mentation is to apply the free propagator Go and to obtain \ ^tpi). Finally (\E'f, | is acted 
on obtaining the transition amplitude N'^^. This method we denote by "method '^nn+znf-^ 

In the absence of the 3N force Eq. (34) simplifies to 

f\(Pi)=tP\ (Pi) + tPGof I (Pi) (35) 

and the corresponding result will be denoted as "3 at at." In our numerical implementation 
we always calculate the matrix elements of the current operator in the frame in which 
the photon momentum is parallel to 2;-axis. To obtain the transition amplitude for each 
final angle between the outgoing photon and the beam direction we have to adjust the 
corresponding initial deuteron or proton beam direction. This demands the repetition of the 
iteration of Eq. (34) in each case. One can avoid this by calculating the matrix elements of 
the current operator in a general frame. 

Of course, the methods '^NN+mr and ?>nn by construction can be used only for the 
Nd-radiative capture and the two-body photodisintegration. 
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III. RESULTS 



We would like to present the quality of our methods for pd-capturc for six exemplifying 
observables: the nucleon and deuteron vector analyzing powers Ay{N) and the 
differential cross sections and the tensor analyzing powers Ayr, Axx and T21 = ^Axz- 

Predictions based on the AV18 interaction for these observables for pd capture are pre- 
sented in Fig. 1. The incoming deuteron laboratory energy is = 300 MeV. This is 
equivalent to a photon laboratory energy pal06 MeV in the photodisintegration process. 
The iteration of the Faddeev-hke equations uses a lot of computer power. Thus in the case 
of method Sjviv which in our numerical implementation, as mentioned above, demands a 
separate solution of Eq.( 35) for each angle of the outgoing photon we present only a few 
points (denoted by crosses). We sec that the methods Inn and 3nn agree nicely in all cases, 
except for the tensor analyzing power Axx at small angles. This agreement is obtained 
taking into account in both methods the partial waves with total angular momenta in the 
two-body subsystem up to j' = 3. The method 2jvAr demands much more partial waves (up 
to j — 5) and still there are more angular ranges where the predictions of method 2^^ differ 
from those of methods Inn and 3nn- As mentioned above this is due to the lack of the 
t-matrix in the driving term. 

A nice agreement is obtained when comparing results of methods Inn+snf, '^nn+snf 
and Snn+3NF- This is shown in Fig. 2 for the same deuteron energy, E(i=300 MeV, and for 
a much lower one, E(i=17.5 MeV, in Fig. 3. The cross section is especially insensitive to 
the used method. In addition to predictions based on the AVIS and the Urbana IX forces 
in Figs. 2 and 3 we present also results based on the AVIS interaction alone. As we see, 
except for the cross section, the three nucleon force effects are very small. Comparison to 
the Sagara data [18] at E^^IT.S MeV shows a reasonable agreement for the tensor analyzing 
powers Axx and Ayr and the known disagreement for vector analyzing power Ay(d) [5]. 
Comparing our predictions to the Pickar data [19] at Ed=300 MeV we see that for the proton 
analyzing power Ay{N) we agree with the data at the lower and middle angles, while at 
the higher angles all theoretical predictions are above the data. In the case of the cross 
section differences are much smaller, however the theoretical predictions seem to be flatter 
than data. It is also very apparent that the inclusion of the Urbana IX 3NF improves the 
description of the data. Only two experimental points at very small and at very large angles 
are in better agreement with the pure NN force prediction. 

As is well known, for low energies the 3NF contributes mainly in the 3N bound state, 
while for higher energies the 3NF effects are seen also in the continuum. This is the reason 
why we present results also for a relatively high energy (Ej=300 MeV). We also note that 
the agreement between all three methods is slightly better for lower energy, especially for 
the deuteron tensor analyzing powers (e.g. Axx)- 

The next two figures show dependences of our results on the number of partial waves 
for the example of method 1nn+3NF- In Fig- 4 we show the convergence of the predictions 
in the total angular momentum Jmax of the three-body system. While using only channels 
with Jmax = I is absolutely insufficient, using channels with Jmax = I very close to the 
final prediction with Jmax = y- Allowing the transitions to higher total angular momenta 
states (up to Jmax = y) '^^^^ change the predictions in a perceptible manner. A similar 
picture appears for the convergence in the number of partial waves used in the two-body 



S 



subsystem (sec Fig. 5). While using only channels with total angular momenta in the 
two-body subsystem up to j = 2 is far from the predictions with j = 4, there is only a 
small difference between the predictions with maximal values of the two-body total angular 
momentum j — 3 and j — 4. However, the method 2]\f]\f (not shown in Fig. 5) requires at 

least jmax — 5. 

Finally, we would like to compare predictions for the three-body photodisintegration of 
■^Hc. In Figs. 6-9 we present examples of exclusive differential cross sections for different 
kinematical configurations at E^=100 MeV, given as a function of the S-curve arc-length. 
The two protons are measured under different polar and azimuthal angles ©, In all 
cases we see that there is an excellent agreement between the predictions based on methods 
IjvTV+sAfF and 27VAr+3iVF- They are represented by solid and thick dotted curves, respectively. 
In both cases we show results with channels up to jmax < 3 and Jmax ^ As we checked 
for a large number of configurations (about 300 000) the differences remain for all cases 
below 1%. In the case of the predictions based on NN interaction only, the differences are a 
bit bigger - for the majority of the configurations they are below 5%. In that case channels 
with Jmax < 4 were used for the method Iattv (dotted line) and with jmax < 5 for the method 
27V AT (dashed line). Again, the reason are the different forms of the driving term in the two 
methods. 



IV. SUMMARY 

Different formulations of Faddeev-like equations for pd-capture (equivalent to two-body 
photodisintegration of the 3N bound state) and for three-body photodisintegration of the 
3N bound state have been investigated. This is important to guarantee reliable and well 
converged theoretical results for modern NN and 3N forces. Such tools allow an unambiguous 
test of the dynamics when compared to data. This study is especially timely since a new 
approach to nuclear forces and currents based on effective field theory constrained by chiral 
symmetry is under vivid development [20] - [21]. As we demonstrated all methods show a 
good agreement. Also the computer time needed is very similar. Therefore, none of them 
is preferable, but all are equally useful. The comparison, however, was a very important 
internal benchmark of the methods. Only the method 2jv7v is much slower, since there 
more two-body partial waves have to be used. For that reason we have not checked another 
fourth possible formulation with the same driving term as in the method 2^^^, but also three 
nucleon interaction in the integral kernel [7] . The reason for the slow convergence is that the 
permutation operator from the integral kernel acts onto the permutation operator from the 
driving term and one has to use a much bigger number of partial waves to obtain converged 
results. 

We can conclude that now different methods are available, which can be considered to 
be very reliable and the achieved results document the numerical accuracy up to the order 
of a few percent slightly dependent on the observable. At the low energies the accuracy is 
much better. In addition this can be improved if smaller experimental errors in the future 
will require that. 
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FIG. 1. The comparison of the three methods using the AV18 interaction at the deuteron 
laboratory energy £^=300 MeV. The sohd hne represent predictions based on method Inn, the 
dashed one on method 2^^ (here jmax=5) and the crosses are for method SjvAr (see text). 
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FIG. 2. The comparison of the three methods using the AV18+UrbanaIX interaction at the 
deuteron laboratory energy Ed=300 MeV, jmax=3, Jmax = The soUd hne represent predic- 
tions based on method 1nn+3NF, the dotted one on method 2nn+3NF and the crosses are for 
method 3nn+3NF- The dashed hne represents predictions of method Inn, based only on the 
AV18 interaction. Data (o) are from [19]. 
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FIG. 3. The comparison of the three methods using the AVlS+UrbanalX interaction at the 
deuteron laboratory energy Ed=17.5 MeV, jmax=3, Jmax = Curves and crosses as in Fig. 2. 
Data (o) are from [18]. 



14 




FIG. 4. The convergence in the number of three-body partial waves for the method 1nn+3NF 
using the AV18-|-TJrbanaIX interaction. E(^=300 MeV, jmax— 3 fixed, Jmax — o 

(dash-dotted), I 

(dotted),^ (dashed), ^ (solid line). Dashed and solid lines are indistinguishable. 



15 




I I I I I I I I I I I I I I I I I 

60 120 180 60 120 



0^ [deg] [deg] 



FIG. 5. The convergence in the number of two-body partial waves for method 1nn+3NF using 
the AV18+UrbanaIX interaction. Ed=300 MeV, Jmax=^ fixed, jmax = 2 (dashed), 3 (dotted) and 
4 (sohd hne). 
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FIG. 6. The differential cross section for three-body photodisintegration at E-y=100 MeV at 
the two proton angles: Gi = 10°, $i = 0°, 62 = 10°, $2 = 0°. The dotted, dashed, solid and thick 
dotted curves represent methods Inn, ^nn, '^nn+3NF, ^nn+3NF, respectively. 
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FIG. 7. The differential cross section for three-body photodisintegration at E^=100 MeV at 
protons angles: Gi = 90°, $1 = 0°, G2 = 90°, ^2 = 90°. Curves as in Fig. 6. 
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